On the critical temperature and the energy gap in dense SiH 4 (H 2 )2 at 250 GPa 
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The critical temperature (7b) and the energy gap (2 A (7 1 )) for the superconductor SiH4(H2)2 at 
250 GPa have been calculated. The wide range of the Coulomb pseudopotential's values has been 
considered: fj,* £ (0.1,0.3). It has been stated that Tc decreases together with the increase of [i* 
from 129.83 K to 81.40 K. The low-temperature energy gap (T ~ K) decreases together with the 
increase of the Coulomb pseudopotential from 50.96 meV to 30.12 meV. The high values of 2 A (0) 
mean that the dimensionless ratio Ra = 2A (0) /fcsTb significantly exceeds the value predicted by 
the classical BCS theory. In the considered case: _Ra £ (4.55, 4.29). Due to the unusual dependence 
of the critical temperature and the energy gap on fi* , the analytical expressions for 7b (/x*) and 
A (/x*) have been given. 
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The high-pressure superconducting state, induced by 
the electron-phonon interaction, may have the very high 
value of the critical temperature. Due to the small mass 
of the nuclei forming the crystal lattice and the lack of 
the external electron shells, the most interesting proper- 
ties should have been revealed by the high-pressure su- 
perconducting phase in the hydrogen 0], [H, 

The ab initio calculations suggest that the metalliza- 
tion of the hydrogen takes place for the pressure p ~ 400 
GPa [f|. In the considered case, the hydrogen exists in 
the molecular phase. It has been stated that the criti- 
cal temperature in the molecular hydrogen reaches high 
values (~ 240 K for p = 450 GPa); the other thermo- 
dynamic parameters are significantly different from the 
predictions of the classical BCS theory @, Q, @- Above 
the pressure equal to 500 GPa follows the dissociation of 
the hydrogen's molecular phase into the atomic phase Q , 
0; 01- The superconducting state in the atomic 
phase has been analyzed up to the value of 3.5 TPa [|[. 
The obtained results suggest the extremely high critical 
temperature (~ 600 K for p = 2 TPa) 0. In the con- 
sidered case, the ratio of the low-tcmpcrature energy gap 
and the critical temperature equals more than 6, so it is 
comparable with the values measured in the cuprates Q , 

Because of the high metallization's pressure, the exper- 
imental analysis of the hydrogen's superconducting state 
is unavailable. For this reason, currently examined is the 
existence of the high-temperature superconducting state 
in the compounds of IV group: CH4, GeH4, SnH4, and 



SiH 4 [13(, [I4j, [15J, M, [17|, UJ, 1190- 

In the case of SiH4, the metallization occurs at the 
pressure 50 — 60 GPa (T^, The low-temperature 

superconducting state (Tc — 17 K) has been observed for 
the pressure at 96 GPa and 120 GPa [l9[. It is possible 
that the higher values of the critical temperature would 
be measured in the range of the higher pressures or in 
the compounds of the silicon like: Si2H@ or SiH| (H2) 2 

In the case of SiHj (H2) 2 , it has been stated that the 
critical temperature can take the value of 107 K for the 
pressure at 250 GPa [2l[. We note that Tc has been 
estimated on the basis of the McMillan formula [22| . 

From the microscopic point of view, the growth of 
the critical temperature in SiH4 (H2) 2 is related to the 
strong interaction between the additional molecules H 2 
and SiH4, which leads to the significant increase of the 
electron-phonon coupling constant. 

In the present study, wc have determined the exact 
dependence of the critical temperature and the energy 
gap on the Coulomb pseudopotential for S1H4 (H2) 2 com- 
pound (p = 250 GPa). Due to the very high value of the 
electron-phonon coupling constant (A = 1.61), the cal- 
culations have been carried out by strictly solving the 
Eliashberg equations [23 1 . 

The Eliashberg equations on the imaginary axis as- 
sume the following form: 
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The quantity </>„ = 4>(iu} n ) represents the order pa- 
rameter function; Z n = Z (iu) n ) denotes the wave func- 
tion renormalization factor; n-th Matsubara frequency is 
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defined by the expression: oj n = (n/ f3) (2n — 1), where 
P = {^bT) , and ks is the Boltzmann constant. In 
the Eliashbcrg formalism, the order parameter is given 
by the formula: A = <fi/Z. 

The symbol A (z) denotes the pairing kernel for the 
electron-phonon interaction: 



cut-off frequency; oj c = 3f2 max . 

The Eliashbcrg equations have been solved for 2201 
Matsubara frequencies (M = 1100). The numerical 
methods used in the paper have been presented in: (24|, 
HH, HH, and 27 1 . The convergence of the solutions has 
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where a 2 F (ft) is the Eliashbcrg function. For SiH4 (H.2) 2 
under the pressure at 250 GPa, the Eliashberg function 
has been determined in the paper plj . The maximum 
phonon frequency is equal to 375.98 meV. 

The depairing Coulomb interaction is parameterized 
with the help of the Coulomb pseudopotential /Lt*. The 
symbol 8 denotes the Heaviside function and oj c is the 



been obtained for T > T = 11.6 K. 

Using the Eliashbcrg equations on the imaginary axis 
one can precisely calculate the value of the critical tem- 
perature in the dependence on the Coulomb pseudpoten- 
tial. In the study, it has been assumed fi* G (0.1,0.3); 
due to the fact that the exact value of the Coulomb pseu- 
dopotential is very difficult to calculate with the use of 
the ab initio methods. 

In order to estimate the physical value of the order 
parameter, the Eliashberg equations should be solved in 
the mixed representation (the equations are defined si- 
multaneously on the imaginary and real axis) (28{: 
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The symbol N (oj) and / (oj) denotes the Bose-Einstein 
function and the Fermi-Dirac function, respectively. The 
equations and (O have been solved with the help of 
the numerical methods used in the papers: [29j], [lo[ and 

mi. 

The solutions of the Eliashbcrg equations on the imag- 
inary axis have been presented in Fig.Q] In particular, in 
Fig.Q] (A)-(C) there is plotted the dependence of the or- 
der parameter on the number m for the selected values of 
the temperature and the Coulomb pseudopotential. On 
the other hand, in Fig.Q] (D)-(F) there is presented the 



influence of the temperature and the Coulomb pseudopo- 
tential on the form of the wave function renormalization 
factor. 

On the basis of the obtained data, it has been found 
that the increase of the temperature causes the significant 
decrease of the maximum value of the order parameter 
(A rn= i); also decreases the half-width of the considered 
function. In the case of the increase of the Coulomb pseu- 
dopotential, the function of the order parameter strongly 
decreases in the range of the lower Matsubara frequencies 
and becomes saturated at the level of the progressively 
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FIG. 1: The dependence of the order parameter (figures (A)-(C)) and the wave function renormalization factor (figures (D)- 
(F)) on the number m for the selected values of the temperature and the Coulomb pseudopotential. The first 200 values of the 
functions A m and Z m have been presented. 



smaller values. 

Analyzing the results obtained for the wave function 
renormalization factor, the very weak dependence of the 
function Z m on T and fj* can be easily observed. From 
the physical point of view, the achieved result indicates 
that the electron effective mass is slightly responsive to 
the changes in the values of the temperature or /i* pa- 
rameter. 
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FIG. 2: The maximum value of the order parameter (figure 
A) and the maximum value of the wave function renormaliza- 
tion factor (figure B) as the function of the temperature for 
the selected values of the Coulomb pseudopotential. 



The temperature dependence of the order parame- 
ter and the wave function renormalization factor can 
be traced in the most convenient way by plotting the 



courses of the functions A m= i (T) and Z m —\ (T). The 
results have been presented in Fig. [2] (A) and (B). In 
the case of the order parameter, the quantity 2A m= i (T) 
with the good approximation reproduces the tempera- 
ture evolution of the energy gap at the Fermi level. Let 
us notice that the shape of the function A m= i (T) can 
be parameterized with the help of the following expres- 
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33.01 meV and 
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In the case of the second solution of the Eliashberg 
equations, Fig. [5] (B) confirms the very weak dependence 
of Z m= i on the temperature and the Coulomb pseudopo- 
tential. Note that for T = Tc, the maximum value of the 
wave function renormalization factor is independent of 
(i*. In the considered case, it can be calculated with the 
help of the simple formula: [Z m= i] T _ T =1 + A = 2.61. 
The identical value has been achieved on the basis of the 
numerical calculations. The above result proves the high 
accuracy of the used numerical methods. 

Below, there is determined the strict dependence of the 
critical temperature on the Coulomb pseudopotential. In 
the considered case, the shape of the function T c (//*) has 
been reproduced numerically for 300 values of fj,* from 
the range 0.1 - 0.3; the condition [A m=1 ] T=Tc = has 
been used. The results have been presented in Fig. [3j In 
addition, there arc also plotted the critical temperature 
obtained using the classical McMillan formula and the 
Allen-Dynes expression [22j], (32|. It is easy to see that 
the critical temperature calculated in the analytical way 
is underestimated, especially in the range of the higher 
values of /i* . However, the Allen-Dynes expression much 
better predicts Tc than McMillan formula. 
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FIG. 3: The critical temperature as the function of the 
Coulomb pseudopotential. 



In view of the difficulties in precise estimating the value 
of the critical temperature by the classical expressions, 
the new formula for Tc has been given. In particular, 
we have used the method of least squares and 300 exact 
values of the function Tc (n*)- The obtained result takes 
the form: 



ksTc = /i/2:j-^exp 



1.125 (1 + A) 

X-fj* 



(6) 



where fi and / 2 denote the correction functions [32J 
and / 2 = 1 + "T2 i a 2 — ■ The quan 
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tities Ai and A 2 have the form: Ai = 2 — 0.14/i* and 
A 2 = (0.27+ 10//*) ( v / W2/win)- The parameters wi n and 
y/U2 are equal to 71.76 meV and 122.67 mcV, respec- 
tively. 

On the basis of Fig. [31 it is easy to see that the formula 
(|()]) exactly reproduces the numerical values of the critical 
temperature. 

The solutions of the Eliashberg equations on the imag- 
inary axis are used as the input parameters to the Eliash- 
berg equations in the mixed representation. Those equa- 
tions allow to continue analytically the order parameter 
and the wave function renormalization on the real axis: 
A m -> A (oj) and Z m ->• Z (cj). 

In Fig. HI the results obtained for the order parameter 
have been presented. In particular, in Fig. [4] (A)-(C) we 
have plotted the courses of the real part of the function 
A (ui) for the selected values of the temperature and the 
Coulomb pseudopotential; ui G (0,w c ). It has been stated 
that in the range of the frequencies from to f2 max the 
function Re [A (w)] has the complicated shape, which can 
be divided into the number of the local maxima and min- 
ima (especially in the range of the low temperatures) . It 
should be emphasized that for ui G (0, fl ma x) the com- 
plex behavior of the real part of the order parameter 
is induced by the complicated shape of the Eliashberg 
function. Analyzing Fig. 2] (A)-(C), it can be addition- 



ally noticed that the increase of the temperature smooths 
the course of the function Re [A (to)}. In the range of the 
higher frequencies (ui G (n max , 0J C )), the real part of the 
order parameter is subjected to the saturation. 

The imaginary part of the order parameter on the real 
axis is plotted in Fig.@](D)-(F). There are clear similari- 
ties between the behavior of the function Im [A (w)] and 
Re [A (u>)]. In particular, in the range of the frequencies 
from to f^max, the function Im [A (w)] is characterized 
by the complicated course, closely correlated with the 
shape of the Eliashberg function. For the higher frequen- 
cies it also becomes saturated. 

In order to perform the thorough analysis, the func- 
tion A (w) should be plotted on the complex plane (Fig. 
[5](A)-(C)). After making the corresponding transforma- 
tions, the characteristic spiral forms have been obtained. 
On the basis of the presented results, the complicated 
structure of the order parameter for the frequencies from 
to f2 max can be seen very clearly. 

The obtained results allow to identify the values of the 
frequency for which the effective potential of the electron- 
electron interaction is attractive (Re [A (a;)] > 0) [HI]. In 
the considered case these may be two compartments: the 
first one extends from to uj p < f2 max and exists up to 
the critical temperature and for G (0.1, 0.3). It should 
be strongly emphasized that the increase of the temper- 
ature and the Coulomb pseudopotential causes the clear 
decrease of w p . The second compartment of the frequency 
forms in the range of the lower values of T and fj,* and is 
related to the right hand part of some loops. 

With the explicit form of the function A (ui) , one can 
calculate the exact value of the energy gap at the Fermi 
level. To do this, it is necessary to use the equation: 



A (T) = Re [A (w = A (T))] . 



(7) 



On the basis of the low-temperature solutions (T = Tq) 
and /i* G (0.1,0.3) it has been stated, that 2A (0) G 
(50.96,30.12) meV. Thus, the dimcnsionless ratio Ra = 
2A (0) /ksTc assumes the high values from 4.55 to 4.29. 
Let us notice that the classical BCS theory predicts 
[-RaJbcs = 3-53 [HI. From the physical point of view, the 
achieved result indicates that the thermodynamic prop- 
erties of the superconducting state in SiH4 (H 2 ) 2 very sig- 
nificantly differ from the properties of the canonical su- 
perconducting state (the BCS theory). 

In order to determine the values of Ra ratio it was nec- 
essary to conduct the complicated and time consuming 
calculations. For this reason, we present the analytical 
formula for Ra , which allows to reconstruct the accurate 
numerical results: 
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0.3447. 
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FIG. 4: The real part of the order parameter (figures (A)-(C)) and the imaginary part of the order parameter (figures (D)- 
(F)) on the real axis for the selected values of the temperature and the Coulomb pseudopotential. The blue and orange lines 
represent the solutions for uj £ (0, f2 ma x), whereas the grey lines correspond to the solutions for uj £ (Sl max ,o; c ). 
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FIG. 5: (A)-(C) The order parameter on the complex plane for the selected values of the temperature and the Coulomb 
pseudopotential. The lines with symbols represent the solutions for uj 6 (0, fi ma x), whereas the regular lines correspond to the 
solutions for uj € (f2 max ,^ c ). 



The expression (|8|) has been constructed basing on 300 
exact values of the function Ra (fi*) in the range of /j* 
from 0.1 to 0.3. 

The second solution of the Eliashbcrg equations 
(Z (oj)) serves as the base to calculate the exact value 
of the electron effective mass (m* ) . Performing the cor- 
respondent calculations, it has been found that the elec- 
tron effective mass weakly depends on T and /i*. The 
quantity m* takes its maximum value at the critical 



temperature: m* = 2.94m e , where the symbol m e de- 
notes the electron band mass. The achieved value has 
been calculated on the basis of the following expression: 
m* e /m e = Rc [Z (uj = 0, T = T c )]. 

In the last step let us notice that the value of the elec- 
tron effective mass accurately reproduces the value of the 
expression 1.131 (1 + A), that is included in the critical 
temperature formula (Eq. ©). 

In summary: we have precisely calculated the value of 
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the critical temperature and the energy gap at the Fermi 
level for SiH4(H 2 )2 compound under the influence of the 
pressure at 250 GPa. We have examined the wide range 
of the Coulomb pseudopotential's values; //* € (0.1,0.3). 
It has been stated that Tc can undergo the change from 
129.83 K to 81.40 K. The obtained result means that 
even for the high value of /z* the critical temperature is 
very high. It should be underlined that the exact values 
of the critical temperature cannot be calculated with the 
help of the classical formulas. 

In the considered range of the Coulomb pseudopoten- 
tial's values, the low-temperature energy gap changes 
from 50.96 meV to 30.12 meV. The high values of 2A (0) 
cause the dimensionless parameter Ra to greatly exceed 
the universal value 3.53 predicted by the classical BCS 
theory. From the physical point of view, the discussed 
result indicates that the thermodynamic properties of 
the superconducting state in SHI4 (112)2 significantly dif- 
fer from the properties of the BCS superconducting state. 

In the presented paper, the analytical expressions for 



Tc (/i*) and Ra (a**) have been given. In the future, they 
will allow to obtain the physical Coulomb pseudopoten- 
tial's value for SHI4 (112)2 compound in the experimental 
way. 
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